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Abstract 

We studied the performance of the multicanonical algorithm by simulating the 
van Hemmen spin glass model and reproduced the exact results for this mean field 
model. Physical quantities such as energy density, specific heat, susceptibility and 
order parameters are evaluated at all temperatures. 

One of the recent and promising developments in computational physics is the use of 
the multicanonical ensemble]]]] for numerical simulations. It was originally developed for 
systems with first order phase transitions to avoid the supercritical slowing down || ^ || . 
Another possible target of the method are systems with conflicting constraints. For low 
temperatures these systems split into many thermodynamic states, separated by high 
tunneling barriers. Low temperature canonical simulations tend to get trapped in one of 
those states. Multicanonical simulations overcome these barriers by connecting back to 
high temperature states. First promising studies exist for spin glasses || ||, 0], random 
ising model || and proteins ||. Here, we test the performance of the multicanonical 
algorithm against exact results for this kind of systems, namely the van Hemmen spin 



glass model |10 



In this paper, due to limited space, we will not review the multicanonical ensemble. 
Interested readers should refer to the references [P], Q [□]]■ 



The van Hemmen spin glass model [TQ] is denned by the hamiltonian: 



# = -^E s(i)s(j) - £ jyswsv) - h £ s®, (i) 

describing N Ising spins interacting with an external magnetic field h and with each other 
in pairs Via Jo a direct ferromagnetic coupling is incorporated. The contain the 

randomness, 

J ij = JfitiVj + tjVil, (2) 
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where the & and rjj are independent and evenly distributed random variables with mean 
zero and variance one. This distribution of Jjj 's is shown |T2| to model the RKKY [|13[] 



interaction in a real metallic spin glass: symmetric and highly peaked at = 0. The Jjj 
contain 2N independent random variables and describe therefore a random-site problem, 
not a random-bond problem like the SK model. 
The model has three order parameters, 

N N N 

m N = N- 1 Y,S(t), q 1N = N- 1 J2^S(t), q 2N = iV" 1 ^^(i). (3) 

i=l i=l i=l 

Without a ferromagnetic interaction, i.e., for Jo = the magnetization vanishes m = 
0, and the order parameters qiN and q2N are combined to give a more relevant order 
parameter Q: 

Q = N^ <-{q 1N + q 2N )>, (4) 

from which the thermodynamical quantities can be obtained. This model is exactly 
solvable without replicas and its main features are consistent with that of a spin glass 
model with randomness and frustration, except for the metastable state structure. The 



system actually picks out a Mattis state [14] and lacks therefore a great multiplicity of 



metastable states which is considered integral to a true spin glass. We like to mention 



that the van Hemmen model is also very closely related to the Hopfield model [15] which 
is widely used to model neural networks. 

We performed multicanonical simulations of the van Hemmen model on cluster of RISC 
workstations. Independent gaussian distributions for £j and r\i with mean value zero and 
variance one were created. Simulations with up to iV = 1000 spins were easily carried out. 
Thermodynamical averages were evaluated over typically two million iterations, following 
2 x 10 5 iterations of thermalization runs. With this extended statistic and the new method 
we could go far beyond what was done in earlier work[[L6|]. Let us first concentrate on 



Figure 1: Spin glass order parameter distribution. 



the pure spin glass case and set the ferromagnetic coupling J = 0. Fig.l shows the 
distribution of the order parameter Q at all temperatures for N = 1000 spins . Tf/J=l 



Figure 2: Specific heat vs. temperature for several lattice sizes. 



is the bifurcation point below which temperature the nonzero ordering sets in and reaches 
its maximum value \jypK at T = 0. There is no many- valley structure like the one 
observed in the multicanonical simulation of Edwards- Anderson model 0. The internal 
energy assumes the values at T = Tf and — 1/n at T = 0. Fig.2 displays the specific heat 
vs. temperature for several lattice sizes. The specific heat is linear at low temperatures, 
peaks at 0.8 < T < 1.0 and vanishes for Tf < T. The linear behavior of the specific 
heat indicates the existence of many low lying excited states at low temperatures. Both 
our values for the internal energy and specific heat reproduce the exact values with high 
precision. 

Next we included the ferromagnetic coupling in the model, with zero external field. 
The spin glass to ferromagnet transition is supposed to take place in the region J ~ J- 
We observed that the magnetization jumps and the order respectively vanishes for Jq 
approaching J. While at Jo/ J = 0.6 the distribution of the order parameter was the same 
as depicted in Fig.l, for Jo/ J = 0.9 it assumed the value Q = 0.34 at about T/J ~ 2 
and stayed constant all the way down to T = 0. The system stays in the metastable 
spin glass phase, and does not spontaneously jump to the ferromagnetic phase. Such 
behavior was observed elsewhere ]I7J. In Fig. 3 we show the zero field susceptibility Xo(T) 
vs. temperature for several values of the ferromagnetic coupling J . For T < Tf the 
susceptibility develops a plateau, a feature indicating an equilibrium phenomena, and also 
shared by the Sherrington-Kirkpatrick model. The value of the susceptibility at T = Tf 
agrees well with the analytic solution Xo(Tf) = l/( J — J ). For Tf < T, we observe a pure 
Curie- Weiss behavior Xo(^) = VC^ ~~ ^7)- Entering the region of the metastable spin 
glass phase mentioned above, around J / J = 0.6, the plateau starts getting distorted. 

The multicanonical simulation of the van Hemmen model is quite feasible; the phases 
are clearly observed, main features like linear rising specific heat, the plateau in zero- field 
susceptibility are easily obtained. Resulting numerical values agree remarkably well with 
the predicted ones. It is important to notice that these thermodynamic quantities were 
obtained for all temperatures from one simulation. In our simulations it is straightforward 
to probe the ground states for systems like the present one which includes randomness 
and frustration and a continuous energy spectrum. The multicanonical algorithm reveals 
itself to be a reliable instrument for the simulation of systems with conflicting constraints 
and a useful tool to study NP-complete systems. 



Figure 3: Susceptibility vs. temperature for several values of the ferromagnetic coupling. 
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